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The advent of high-temperature superconductors has re-invigorated the study of lay- 
ered superconductivity, 1-3 wherein adjacent layers are Josephson coupled. An issue that 
remains open is the question of whether or not layer decoupling occurs in the presence 
of a parallel magnetic field. 4-9 The standard phenomenological model used to study this 
problem theoretically is given by the Lawrence-Doniach (LD) free-energy functional, which 
in the absence of fluctuations of the magnetic field reads 9 
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Here, 6>(r, /) denotes the phase of the superconducting order parameter, which lives on 
iV equally spaced square lattices with respective coordinates / and r = (x,y), while the 
parallel magnetic field = ($o/27r<i)&|| directed along the y axis is related to the vector 
potential above by A z = — b\\x, where d represents the spacing between layers. Note that 
V = (A x , A y ), with A^6(r) = 6(r + p) — 6(r). Also, J\\ is a measure of the in-plane phase 
rigidity, while 7* is directly related to the mass anisotropy parameter 7 = (mj_/m||) 1//2 of 
the superfluid motion within (my) and between (mi) planes. For the bulk case, — > 00, 
where magnetic screening effects must be included, it was first claimed by Efetov 4 that 
layers decouple in the mixed phase for parallel fields B\\ above the characteristic scale 
Bl ~ $o/7^ 2 - Recent work, however, finds no evidence for a true phase transition as a 
function of field. 6-9 

The author has recently studied the relatively simpler problem of layered thin films 
of extreme type-II superconductors (Az, — > 00) in parallel magnetic field. 9 In this limit, 
magnetic screening effects are negligible and the problem reduces to the study of N weakly 
coupled square-lattice XY models in the presence of uniform frustration; i.e., the thermo- 
dynamics is determined by the energy functional Exy = — J2 r n J^ cos [^n < P( r ) — A^r)], 
where J x = J\\ = J y and J z = J\\/y 2 , with 7' = jd/a, and where = (0, b±x, — b\\x). 
Here a denotes the square lattice constant. For thin enough films, Nd 'j'd Az, in 
perpendicular magnetic fields, 10 B± = ($o/27ra)6j_ ^> $o/7 2 ^ 2 5 the thermodynamics of 
this model factorizes into parallel and perpendicular pieces that correspond respectively to 
Josephson vortices in between layers and to 2D perpendicular vortices within each layer. 9 
The latter is corroborated by Monte-Carlo simulation. 11 Physically, this factorization is 



due to the fact that the modified Josephson penetration length j'd sets the minimum 
perpendicular size for vortex loops that traverse many layers. 12 The present factorization 
prevails in the presence of both parallel and perpendicular magnetic field as well, with 
the parallel thermodynamics determined by the following partition function for a layered 
Coulomb gas ensemble: 9 ' 13 
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where n z (r, I) is an integer field over the layered structure that describes inter-layer fluxon 
excitations, 14 with the fields at the boundary layers set to n z (f, 0) = = n z (r,N), and 
where = — V -2 is the Greens function for the square lattice. Here, (3\\ = J^/ksT, 
while (3±_ = ,9||/7 /2 . (Note that the assumption that magnetic screening be absent requires 
that the perpendicular component of the magnetic field exceed $oMi-) By extending 
Polyakov's Hubbard- Stratonovich transformation of the neutral Coulomb gas in the plasma 
phase, 15 it can further be shown 9 that this layered Coulomb gas is equivalent to the LD 
model (1) at temperatures below the zero-field decoupling transition temperature, T* = 
47rJ||, in the limit of small fugacity, yo = exp(— 7 /2 /2/3||), where the effective anisotropy 
parameter is given by 7^ = (f3\\/2yo) 1 / 2 . 

In this paper, we will first establish that the above layered Coulomb gas ensemble 
(2) is equivalent to a one-dimensional (ID) fermion analogy consisting of coupled chains. 
The latter is a generalization of the repulsive back-scattering model for spin-1/2 fermions 
introduced by Luther and Emery (LE), 16 where the spin is identified with the layer index. 
The Hamiltonian for this model is divided into two parts, H = H\\ + H±, with 
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and with field operators ^ R (x,l) = J2k elkXa (k, I) and ^(x, /) = Lx ^ 2 J2k ^ kx b(k, Z)| 

for right (R) and left (L) moving fermions. We thus have N Tomonaga-Luttinger chains, 
with adjacent chains coupled via a LE-type repulsive back-scattering interaction (U± > 0). 
The magnetic flux between consecutive layers in the LD model (1) is given by 6y = 
27r(iV;_|_i —Ni)/L x , where iVj denotes the number of spin-less fermions in the Z th chain. The 
above fermion analogy thus completes a triad of equivalent descriptions (1), (2), and (3a,b) 
that generalize the known equivalences between the sine-Gordon model, the 2D Coulomb 
gas, and the massive Thirring/LE models to iV layers. 14-18 By analyzing the double-layer 
case (N = 2) that corresponds to the original spin- 1/2 back-scattering model, 16 we find 
that the equilibrium magnetization for parallel fields shows an anomalous B^ 1 tail in the 
high-field limit at temperatures near T* (see Fig. 1), as opposed to the B^ 3 tail expected 
from Ginzburg-Landau theory. 5 ' 9 This is a result of the entropic presure between neigh- 
boring Josephson vortices, which the fermion analogy correctly accounts for. Nevertheless, 
we continue to obtain at best only a cross-over behavior as a function of parallel magnetic 
field (see Fig. 1), which had been claimed earlier on the basis of a semi-classical analysis. 9 

Equivalence. We now proceed to show that the above ID fermion back-scattering 
model is equivalent to the layered Coulomb gas ensemble (2) by extending the demon- 
stration given by Chui and Lee 17 for the minimal case N = 2. The S matrix, which 
measures the overlap of the unperturbed (U±_ = 0) groundstate $o with the exact one 

is (^ol^o) = (^o l ex P( — iToH±) |3>o)- After making the standard canonical transfor- 
mation to trivialize the U\\ term in (3a) and employing the boson representation for the 
back-scattering interaction, 16 ' 19 H±, a perturbative expansion of the former exponential 
yields 
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where n z (xi, Ti, h) = ±1 are inter-layer fluxon charge distributions 13 such that J2i n z(xi, Ti, h) =| 
0, where A z (x,t,1) = —2(kF,i+i — kF,i)x is a function of the Fermi wavenumbers, kF,i 



corresponding to each chain, /, and where the average (...)o is over the non-interacting 
groundstate. Above, <j>j(x, r, /) = ipj(x, r, / + 1) — ipj(x, r, /), where ipj(x, r, /) is the time 
evolution of the local operator ipj(x, I) = \im a ^o2TcL~ 1 J2 k /c _1 exp( — \a\k\ — ikx)pj(k, I) 
that results from the boson representation, with the usual particle-hole operators given by 
PR,(k, I) = Yliq °^ (q + Z)a(g, /) and /) = ^ (q + k, l)b(q, I). Notice that n z = 1 

corresponds to the choice of the step-up term for H± in the perturbative expansion, while 
n z = —1 corresponds to the step-down term. Also, we have that tanh2<^> = U\\/2ttvf, as 
a result of the above-mentioned canonical transformation, 16 which yields a renormalized 
Fermi velocity equal to v' F = vf sech2^>. After separating ipj into creation and destruction 
pieces, changing variables to y = v' F r and L y = v' F iTo, and extending the Chui and Lee 
procedure 17 to the present case (4), we obtain the equality ($ |^o) = -^cg along with the 
following identifications: 

b\\=2(k F ,i+i-k F ,i), (5) 
Z^ 1 =47re 2 ^, (6) 
y =(^)-\a/af(\U ± \/v' F ). (7) 

Notably, since Ni = iv~ 1 kF,iL x gives the number of spinless fermions in a given chain, 
we obtain the relationship b\\ = 2tt(Ni + i — N{)/L x for the average magnetic induction in 
between layers I and / + 1 announced in the introduction. Notice that the back-scattering 
term (U±) sets the anisotropy parameter, 7^, in the LD model (1), whereas the trivial 
intra-chain interaction (U\\) sets the temperature. Finally, if we approximate the S matrix 
of the fermion analogy by ($o|^o) = exp{— iTo[Ef(U±) —Ef(0)]}, where Ep(U±) denotes 
the energy of the groundstate then we can identify the free energy G s — G n of the LD 
model in parallel field (1) with the latter via 

(L y /v' F )[E F (U ± ) - E F (0)] = (G a - G n )/k B T. (8) 

Thermodynamic properties of the layered superconductor are in this way directly related 
to ground-state properties of the fermion analogy. 

Double Layer. Consider now the minimal N = 2 case, where the LD model (1) re- 
duces to the sine-Gordon model, the ensemble (2) reduces to the conventional 2D Coulomb 
gas, and where the fermion analogy (3a,b) is simply the original LE back-scattering 
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model for spin-1/2 fermions. In this instance, charge and spin particle- hole operators 
p'jik) = 2- 1 / 2 [p' J (k,l) + Pj(k,2)] and a'^k) = 2- 1 / 2 [p' J (k, 2) - p'^k, 1)] can be defined, 
which separate out of the canonically transformed Hamiltonian fZ| . Furthermore, at the 
special point 2 1 l 2 e^ = 1 identified by LE, the spin part becomes equivalent to the one-body 
Hamiltonian 

fc A; fc 

for spinless fermions, where the spin particle-hole operators are given by cr' R (k) = J2 q a\, +q a q \ 
and cr' L (k) = J2 q b\ +q b q - Above, A a = U±(27va)~ 1 is the spin gap, while p = \{pt>i — Pi) 
is the external field. We therefore have energy eigenvalues Sk = ±(v F k 2 + A 2 ) 1 / 2 for such 
fermions. Notice that the procedure followed here is to first obtain the spinless fermion 
Hamiltonian in the absence of field, and to then add the trivial chemical potential shift 
in the presence of field. Let us begin by computing the lower-critical field of the dou- 
ble layer superconductor, and hence turn off the external field within the fermion model 
(p = 0). The edges of the double layer contribute Dirac Eq.-type bound states at zero 
energy that decay as e ±x / A - 7 , where Xj = v' F /\A a \ is the effective Josephson penetration 
length. 20 Then Eq. (8) implies that the line-tension of a single Josephson vortex is given by 
£|l = kBT\A a \/v' F = 27rJ||/Aj, since |A CT | gives the energy cost of adding one fermion/flux 
quantum to the system, and since ksT = 2ttJ\\ at this special point [see Eq. (6)]. The 
lower-critical field is then simply = 4-7r£||/$o- Last, we note that the identification (7) 
implies that the bare Josephson penetration length corresponding to the LD model (1) is 
equal to 7^a = (ttv ' F /\U±\) 1 / 2 a at the present temperature, while Xj = 2ir(v' F /\U ±\)a. In 
the limit of weak coupling, U± — > 0, the latter renormalized Josephson penetration length 
is then much larger than the corresponding bare scale. In particular, the ratio of the paral- 
lel lower-critical field to the bare one set by the LD functional is equal to 'fia/Xj ~ afarf*, 
which is small if a ~ a. This effect is a result of vortex wandering. 9 ' 22 

To compute the parallel equilibrium magnetization deep inside the mixed phase, B\\ ^> 
H^, we first relate the external field p to the parallel magnetic induction in the absence 
of inter-chain/layer coupling (U± = 0); i.e., Pauli paramagnetism gives bu = 27r(A^ 2 — 
N\)/L x = 2-rrxoPi where xo = 4:(2ttv' f )~ 1 . Since the magnetization is generally given by 
My = — Qjq[(G s — G n )/L x L y d], where H\\ = B\\ — AnM\\ = (^o/2nd)h\\ is the magnetic 
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field, we have by Eq. (8) that 

-4ttM|| = (4n/^ )(2nk B T/v , F )(d^/dh ll )-^-{[E F (U ± ) - E F (0)]/L x }. 

But dE F /dn = — (iV~2 - iVi) = -L x xo(/« 2 - | A a) 1/2 > w here we have used iV 2 - iVi = 
2 1 / 2 k F L x /n in conjunction with 2 1 / 2 / u = £fc F . If we presume that dji/dh\\ above is given 
by dfj,/db\\ in the absence of inter-layer coupling, we obtain 



-47rM|i = 2~ 1/2 H l } 1 




(10) 



for the equilibrium magnetization, 21 where B* = (2 1 / 2 /7r)($ /Aj<i) is the parallel cross- 
over field. [The ratio of this cross-over scale to the bare one B ~ $0/7* a<i set by the LD 
functional is again B*/Bq ~ a/a^'^, and hence small if a ~ a.] This simple functional form 
is plotted in Fig. 1. Notably, (10) predicts an anomalous — AttM^/H^ = 2~ 3 / 2 B*/B\\ tail 
for relatively high fields B\\ > B*. Straight forward minimization of the LD functional (1) 
in parallel magnetic field leads to a much weaker B^ s dependence at high fields, 5 ' 9 with 
a more pronounced cross-over at This indicates that the entropic pressure generated 
by fluctuations of the vortex array is not negligible at temperatures in the critical regime, 
which is defined by the relationship /c B T ~ J\\(T) (see ref. 9). 

What, however, is the nature of the vortex array itself? Clearly, the mixed phase 
corresponds to a gapless Luttinger liquid in terms of the LE model. Voit has shown that 
algebraic spin-density wave (SDW) correlations in space dominate in such case for U± > 
0. 19 Hence, by the equivalence demonstrated above, this implies that the array of Josephson 
vortices displays algebraic long-range order transverse to the applied magnetic field. This 
physics is accurately captured by the SDW meanfield theory for (3a,b) characterized by 
the self-consistency equation \i = L~ 1 ^ k {c^(k, l)b(k, /)) and gap equation A; = U\\xi + 
U±(xi+i + Xi-i)- A standard analysis of this theory 23 for the case N = 2 yields a phase 
diagram identical to that derived via renormalization group arguments; 17 i.e., an SDW 
phase (%2 = ~Xi) exists for (repulsive) U± > U\\. Notice that this implies a leading 
dependence /c_bT* = 47rJ||[l + \{a / 'cr^) 2 ] for the decoupling transition temperature as a 
function of the anisotropy parameter by Eqs. (6) and (7). If an external field > |A;| is 
introduced, however, a unitarity catastrophe occurs, wherein \Ai\ acquires an imaginary 
part. 23 This we take to be a signal of the Luttinger liquid {k F ^ 0) "instability". 



N Layers. Given that the above meanfield theory works quite well for the minimal 
case N = 2, it should be accurate at least with respect to the phase diagram for general N. 
For large N, an SDW phase xi+i = ~Xi is stable for (repulsive) U± > |^||- 23 We therefore 
argue by comparison with the match in the case of N = 2 that an anomalous B^ 1 tail 
appears in the equilibrium parallel magnetization as well at high fields B\\ > And that 
the lower-critical field is given by the double-layer result as well can be shown by observing 
that flux penetration must begin at the edges, each of which reduce to a double layer at 
low fields and for weak inter-layer coupling. The latter can be seen via the Coulomb gas 
analogy (2), where it becomes obvious that the effect of neighboring layers is to produce 
a weak dielectric renormalization in the effective 2D Coulomb gas that resides at the edge 
in the limit of weak inter-layer coupling. 22 

In summary, we have obtained a fermion analogy (3) for extreme type-II layered 
superconductors in parallel magnetic field (1). In the minimal double-layer case, we predict 
an anomalous tail for the parallel magnetization in the limit of high parallel magnetic 
field, for temperatures in the critical regime. A mean-field treatment of the fermion 
analogy for the Allayer case (3), however, suggests that this dependence is generic to 
extreme thin-film geometries. 23 

The author is indebted to E. Rezayi for discussions. He also is grateful for the hos- 
pitality of the Theoretical Division at Los Alamos National Laboratory, where this work 
originated. This work was supported in part by National Science Foundation grant DMR- 
9322427 and the Spanish Ministry for Science and Education. 
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Figure Captions 



1. Shown is the equilibrium parallel magnetization characteristic of a double-layer su- 
perconductor deep inside the mixed phase (Bn ^> H^) at temperatures in the critical 
regime. Notice that the pronounced B^ 1 tail resulting from entropic pressure practi- 
cally destroys the decoupling cross-over at £?* (compare with ref. 9). 
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